for k=1:100
    for i=1:100
        for j=1:100
            A(i,j) = i*cos(k/100*pi)+j*sin(k/100*pi);
        end
    end
    E=expm(A);
    imagesc(horzcat(A,E));
    colormap(flipud(gray));
    axis equal;
    M1(k) = getframe;
end

movie(M1,3,30);